Quantum phase diffusions of a spinor condensate 
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We discuss the quantum phases and their diffusion in a spinor- 1 atomic Bose- Einstein condensate. 
For ferromagnetic interactions, we obtain the exact ground state distribution of the phase fluctu- 
ations corresponding to the total atom number (N), the magnetization (Ai), and the alignment 
(or hypercharge) (Y) of the system. The mean field ground state is shown to be stable against 
these fluctuations, which dynamically recover the two continuous symmetries associated with the 
conservation of N and M as in current experiments. 
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Since the observation of Bose-Einstein condensation of 
trapped atomic clouds Q, @ , the coherence properties of 
the condensate has become the focus of many theoretical 
studies US 110,111113. Within the mean-field theory, 
it is commonly assumed that the condensate can be de- 
scribed by a U(l ) symmetry breaking field equiv- 
alent to a coherent state assumption of the ground state, 
(see Refs. for discussions of [/(l)-symmetric ap- 

proaches). Although quite successful in providing the- 
oretical understanding to many experimental observa- 
tions, such a coherent state assumption is not necessarily 
consistent with real experimental situations, where the 
fluctuations of the atom numbers are difficult to control 
[l2l 1113 ] . A coherent state leads to a Poisson distribu- 
tion of atoms. For a ground state with average of N 
atoms, the associated number fluctuations are of the or- 
der VN. As was initially pointed out in Refs. H Q , this 
number fluctuation of a coherent state condensate leads 
to the "diffusion" (or spreading) of its initial phase, fn 
a scalar condensate, this diffusion, a dynamic attempt 
to restore the /7(f) symmetry of the interacting atomic 
system, can be studied in terms of a zero mode, or the 
Goldstone mode of the condensate 0. More physically 
meaningful discussions in terms of the relative phase of 
two condensates were studied soon afterwards [lEHlH]- 
Experimentally, starting with the remarkable direct ob- 
servation of the first order coherence in an interference 
experiment ^4|, direct correlations between number and 
phase fluctuations were observed with a condensate in 
a periodic potential [l^, and more recently, in the re- 
markable Mott insulating state obtained by loading 
a superfluid condensate into an optical lattice [Tj. 

The emergence of spinor- f condensates [Tg, 13 (°f 
atoms with hyperfme quantum number F = 1 ) has cre- 
ated new opportunities to understand quantum coher- 
ence and the associated number / phase dynamics in a 
three component condensate pol l2l|. In this paper, we 
investigate the quantum phase dynamics of a spinor- f 
condensate due to atom number fluctuations. We will 
focus on ferromagnetic interactions, when the conden- 
sate wave functions for the three spin components share 
the same spatial mode |23; in this case, the phase fluctu- 
ations translate into fluctuations of the direction of the 
macroscopic condensate spin. 

We consider a system of N spin- 1 bosonic atoms in- 
teracting via only s-wave scattering [2! 0, HI ■ A 



weak magnetic field B (as always exists in an experiment) 
fixes the quantization z-axis such that quadratic Zeeman 
effect can be neglected. In the second-quantized form, 
the Hamiltonian is I2HI2HI 



H = 



2V72 



' 2M 



V ext (r) - huj L F 2 
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where ipj(r) (j = +,0,—) denotes the annihilation op- 
erator for the j-th component of a spinor-f field. The 
trapping potential V ey ± (r) is assumed harmonic and spin- 
independent. The Larmor precessing frequency is ldl — 
B/ig/H with \xb the magnetic dipole moment for state 
\F = l,Mp = f). The pseudo potential coefficients are 
c = 4vr?i 2 (ao + 2a 2 )/3M and c 2 = 47rS, 2 (a 2 - a )/3M, 
with ao (02) the s-wave scattering length for two spin- 1 
atoms in the combined symmetric channel of total spin 
(2). M is the mass of atom, and F is the spin 1 oper- 
ator [23. Hamiltonian (|TJ is invariant under U(l) gauge 
transformation e %e and SO(3) spin rotations U(a, f3, r) = 
e -iF*<* e -iFv0 e -iF*r ( for 5 — q) 23]. A non-zero B or 

the conservation of magnetization M. = N + — Y_ reduces 
the SO(3) to its subgroup SO(2) generated by e~ tF * a . 

Following the Bogoluibov theory we assume there ex- 
ist three 'large' condensate components <pj around which 
we study the small quantum fluctuations (off-condensate 
excitations) via 



(2) 



Nj is the number of condensed 



where N J dr\4> 3 \ 2 
atoms in component j. At near zero temperatures, higher 
than quadratic terms of 5ipj{r) in the Hamiltonian 
are neglected. The assumed mean field ground state 
\N(j)j(f) breaks both continuous symmetries U(l) and 
SO (2), thus we expect to observe multiple zero energy 
Goldstone modes y|. 

In the first order, we obtain the usual coupled Gross- 
Pitaevskii equation (GPE) for the condensate modes 
4>jif)- The quantum fluctuations Sipj obey the usual 



2 



Bogoluibov-de Gennes equations (BdGEs) 27]. The con- 
densate number fluctuations, can be studied through the 
number fluctuation operators Q 



Pj = df L* (f ) Sijjj (f) + fa (f ) (r) 



df[V>t (f)Vj (r) - A^* (f)fa (r)] . (3) 



Using the GPEs and the BdGEs |23| their equations of 
motion are H, 



-i7LP /2, 



from which we find that both P t( 



(4) 

P + P- and 



P + — P_ are constants of motion, an obvious outcome 
since the Hamiltonian (JIJ commutes with operators of 
both the total number of atoms and magnetization. 
We define the phase operator 



ih / df 



(5) 



with 9j(r) the associated phase mode functions. De- 
note Qtot = Q+ + Qo + the canonical quantiza- 
tion condition [Qtot,Ptot] = iti is satisfied if the con- 



straint J + + Jo 
fdf[9*(f)fa(f) 
densate 



J_ = 1 is enforced, where Jj = 
3.c] . As for a scalar or binary con- 



dQ t 



dt 



NuP, 



(6) 



We note that [Qj,Pfc] = ihS jk Jj, [Sip^f), P k ] = 
8jk4>j(f), and [Sipj(f),Qk] = —ih5jk9j(f). These lead to 
5ipj(f) w e^Pj/Jj + fairjQj/ihJj Hill, which helps 
to define the complete dynamic equations for the number 
and phase fluctuations |2S|. 

Including both the initial and time phases, the ground 
state wave functions can be generally expressed as 
\4>j {f)\e~ l ^~ : ' hWL " >t / h+laj with a common chemical poten- 
tial /x, a Lagrange multiplier enforcing the conservation 
of atom numbers. Minimization of the total energy Eq. 
© leads to a+ + a_-2a = (29|. It was further proven 
in Ref. that the steady state solution takes the same 
spatial mode for each of its three spin components (for 
ferromagnetic interactions), i.e. 



fa(f) = ^/n~ faf)e l 



(7) 



where the real-valued mode function <fi(r) is normalized 
to unity, and governed by an equivalent scalar GPE 



h z v 2 

2M 



+ V cxt (r) + c+N4, 2 {r) 



(r) = M 0(r), (8) 



of a scattering length a 2 (note c+ = cq + C2 oc 02). Define 
rij = Nj/N and m = A4/N as the relative atom numbers 
and relative magnetization, n± = (1 ± m) 2 /4 an d rep = 
(1 — m 2 )/2 in the ferromagnetic ground state |22l l2q. 



Similarly, the phase functions share the same spatial 
mode and can be generally expressed as 

dj{r, t) = v ^0(f)e^- J * Rwi W R+io ' , (9) 

with 9(f) governed by the following equation 

h 2 v 2 



2M 



+ Kxt - H + 3c+Nfa 



9(f) = Nufaf). (10) 



The normalization J d f<f> (f) 9(f) = 1/2 determines the 
Goldstone parameter u, which in the Thomas-Fermi (TF) 
limit is u = [c + /(47rag /3)](a h0 /15a 2 ) 3/5 A r ~ 3/5 , with a ho 
the ground state size (Ti/Muj^o) 1 / 2 of a harmonic trap 
with a frequency Who = (^x^y^z) 1 ^ 3 ■ 

Equations ((HJ and I jlOjl are identical to their counter- 
parts for a scalar condensate 3]. In the ferromagnetic 
ground state, individual atomic spins align along the 
same direction, thus they collide only in the symmetric 
total spin F = 2 channel. To conserve the magnetiza- 
tion with respect to the P-field direction, the condensate 
spin direction is simply tilted at an angle 9 = cos -1 m. 
If this direction is taken as the quantization axis, then 
the spinor condensate behaves essentially as a scalar one 
[Uliilll when the P-field is zero. 

With the mode functions for atom number and phase 
fluctuations, we find the zero mode dynamics can be ex- 
pressed in terms of an associated Hamiltonian. After 
laborious calculation we arrive at |28( 



N 



p'Mp'+q' J Bq' 



(11) 



with p' T = (p' + ,p' Q ,p'_), q' T = (q' + ,q' Q ,q'_), redefined 
as canonically conjugated variables p'j = Pj and q'j = 
Qj/Jj. A & B are two Hermitian and positive-definite 
matrices involving only parameters of the system. Thus 
Pzcro _ 0, the assumed ground state, a mean field sym- 
metry breaking state with coherent condensate ampli- 



tudes vN(f>j(f) is stable. The associated quantum fluc- 
tuations of atom numbers can be studied with the lin- 
earization approximation Eq. (J2J. Matrices A & B 
are simplified with O^e = / df(j) 2 (f)9 2 (f) and = 
J df(/} 4 (f). We also find u = e + Ac+O^e with Ne — 
2 J df9(f) [£ — fi + c + N(f) 2 (f)] 9(f) a non-negative quan- 
tity. 

We note that matrix B is diagonalized by an orthog- 
onal transformation q^ = (q' + + q' + q'_)/V3, qu = 
(q' + - q'_)/V2, and q Y = (q' + - 2q' Q + q'_)/V6 31]. The 
corresponding p?j, pm, and py are respectively fluctua- 
tions of the total number of atoms, the magnetization, 
and the alignment. With these new collective operators, 
the zero mode Hamiltonian becomes 



Hy 



N 



ap 2 N + h v\i + c Py 

+ap N p Y + PpmPy + IPnPm + Wy, ( 12 ) 



where all the coefficients are listed in the Appendix. Re- 
placing pj by —ihd/dqj, Hamiltonian (|12fl leads to the 
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FIG. 1: The N-dependence of fi and 3n5Ay for a 87 Rb 
condensate in a cylindrically symmetric trap with u x — uj y — 
lo± = (2tt)100 (Hz) and uj z = Xlo±. A = 0.1 (solid line), 
A = 1.0 (dashed line), and A = 10.0 (dash-dotted line). 



ground state distribution of fluctuations 



with a width 



: exp 



2h 2 Al 



(13) 



YO , 



A 2 - 



3«o 



4c 2 a 



-C 2 00,£ 



1/2 



and the zero point energy E = Ml/2 = TV[— c 2 00^(£ 



4c 2 0^)] 1/2 - We find n 2 A 



2 



V7/10/3 and 



oc 



-(c 2 /7rag )(a h0 /a 2 ) 3/5 7V 2/5 in the TF limit (and also 
taking the small e — ► 0). Figure shows selected results 
for the TV-dependence of £1 and Ay . While f2 increases 
monotonically with TV, the ground state width A Y0 for 
qy peaks at some intermediate values of atom numbers 
and eventually saturates in the limit of large TV (or strong 
interactions). 

The above ground state distribution of H zcm can eas- 
ily be understood. For a condensate with fixed total 
number of atoms and magnetization, the conservations 
of TV and M require (p 2 N ) — {p 2 M ) — 0, which lead to 
(ill) — (?if) = °°> ^- e -' completely diffused phases; The 
ferromagnetic interaction, nevertheless, prepares a cor- 
related ground state such that both distributions for py 
and qy take a Gaussian form with (p Y ) — l/^A.f'o and 
(q' Y ) = fr 2 Ay /2. Such a distribution in py and qy is 
in fact the minimal uncertainty coherent state consis- 
tent with the symmetry breaking ground state. How- 
ever, it is experimentally difficult to produce such a 
condensate having fixed atom numbers and phase fluc- 
tuations. For any initial state, the dynamic solution 
x(t) = T(t)x(0) as governed by the Hamiltonian Eq. 



CH for x T (t) = [pN(t),PM(t),PY(t),q N (t),q M (t),qY(t)} 
can be used. This solution has a simple structure [28j; 
in addition to oscillating terms of the forms cos fit and 
sinfii, the phase fluctuations of q^ and qu also contain 
diffusion terms proportional to Nt, which indicates the 
linearization approximation of Eq. (J2J is valid only for a 
finite duration. 

Finally, let's consider the diffusion of the direction of 
the condensate spin f(f, t) = J^ij w[ tWij^jir, Be- 
cause all three spin components share the same spa- 
tial mode function 22], individual spins of Bose con- 
densed atoms are parallel for ferromagnetic interactions, 
i.e. they act as a macroscopic magnetic dipole pointing 
along the same direction (independent of the spatial coor- 
dinates): fo(r,t) cx (VT — m 2 cos O, — VT — m 2 sin O, m) 
while undergoing precessing with respect to the z-axis. 
= a + — ao + u>Lt. As the phase dynamics attempts 
to restore the U(l) and SO(2) symmetries of the system, 
their respective initial values become irrelevant. The fluc- 
tuation becomes 5f(t) = J df[f(f,t) — /o(r*,t)] with a 



zero average. Using <r^, a 2 M , and a Y to denote the initial 
variances of TV, M, and Y, we find in spherical coordi- 
nates (r,M) (note sin(9 = Vl-m 2 ), Sf r (t) = Sf r (0) 
and Sfg(t) = Sfe(Q). Both are fixed constants 



[&f r (t)f) = 3a 2 Nl 
[Sfo(t)} 2 ) = 



(14) 



1 — rri 



2 ( 3 ™ V 



N 



•2oi), (15) 



(due to conservations of TV and M.) when respective fluc- 
tuations in TV, M., and Y are uncorrelated. In the az- 
imuthal <f> direction a simple phase diffusion results 



W) = 6U(0) - Rdt, 



(16) 



with 5/^(0) = Vl - to 2 [V2q M (0) + V6mq Y (0)] /2h, 
and a diffusion rate Rd = 2eNSfg(0)/fr, proportional to 
Sfg(0), (which can be explained in terms of the single axis 
twisting of T 2 in the isospin subspace of a spior-1 con- 
densate [3l| ). The TV- and A-dependence of the diffusion 
parameter TVe are shown in Fig. |2Jfor a 87 Rb condensate. 
In the TF limit and applying the results of Ref. [32j , we 
find that e = (Ah 2 /5MR 2 ) ln(i?/1.3a ho ) oc TV" 2 / 5 , gener- 
ally a very small quantity. R = ah (15a 2 /ai 10 ) 1 ^ 5 TV 1 / 5 is 
the TF radius. Typically, the dephasing rate is a fraction 
of uj± for a condensate of ~ 10 5 atoms and A ~ 1, i.e. its 
macroscopic spin direction is lost in a few cycles of trap 
oscillation. Figure|21also indicates that the faster a ferro- 
magnetic condensate diffuses its pointing direction along 
the precession direction, the larger and more tightly con- 
fined a condensate is. 

When B = 0, zero mode dynamics becomes much sim- 
pler if the condensate spin direction is taken as z-axis. 
Although there still exist more than one Goldstone mode 
in this case, the mean field ground state corresponds to 
all atoms in state |+) |2j |3(j . Thus to leading order, 
the only relevant fluctuation is that of the total num- 
ber of atoms (corresponding to the U(l) symmetry), and 
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FIG. 2: The same as in Fig. [I] but for Ne. 



the phase diffusion dynamics becomes the same as in a 
scalar condensate [3j. Our theory as developed in this 
paper does reduce to this simple limit. 

In conclusion, we have studied in detail quantum phase 
diffusions of a spinor-1 condensate with ferromagnetic in- 
teractions. The condensate ground state is very simple, 
all three spin components have the same spatial mode 
function and their associated phase functions are also 
identical. We have constructed the zero mode Hamilto- 
nian for the condensate number and phase fluctuations, 
and solved for the ground state distribution of these fluc- 
tuations when both N and M. are conserved. Further- 
more, we have obtained analytically the dynamic num- 
ber and phase fluctuations relating to both the quantum 
phase diffusion and the initial distribution of these fluc- 
tuations. We have identified a quantum phase diffusion 



coefficient for the pointing direction of the condensate 
spin and recovered its small-time quadratic t-depcndcnt 
spreading. 
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APPENDIX A: COEFFICIENTS IN EQ. (12) 

The coefficients in the zero mode Hamiltonian (112(1 are 



4n 2 



b = 



-(5 + 3m 2 ) + 2O 0e [3c o (l-m 2 ) 2 



+ -c 2 (l-3m 2 ) 



—2 [e(l + m 2 ) - 8c 2 < t >e m 2 



7 

n 



^(l + 3™ 2 )( £ _4c 2 

T\o [e - c 2 <IWl + 3m 2 )] , 

tIq v «j 

-3c 2 n 2 l £>04,/?i 2 . 
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